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ABSTRACT 

The second yntit of th« proposed three year research project has 
now elapsed. The computer automated data acquisition system for atmos* 
pheric emittance, and global solar, downwelled diffuse solar, and direct 
solar irradiances are now fully operational. Hourly-integrated global 
solar and atmospheric emitted radiances have been measured eontinuously 
from February 1981 and hourly-integrated diffuse solar and direct solar 
irradiances have been measured continuously from October 1981. One- 
minute integrated data are available for each of these components from 
February 1982, All Phase II objectives related to the data acquisition 
system have been attained. The results of the correlation of global 
insolation with fractional cloud cover for the first year's data set 
have been completed. A theoretical model has been developed for correlation 
of satellite radiometric data with insolation measurements at the ground. 

A February data set, composed of one-minute integrated global insolation 
and direct solar irradiance, cloud cover fractions, meteorological data 
from nearby weather stations, and GOES East satellite radiometric data, 
has been collected to test the model and develop the cloud dependence for 
the local measurement site. 
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INTRODUCTION 


A solar enargy tneasurement station was established at Hampton Institute, 
February 15, 1981. Routine measurements are being made of global, diffuse, 
and direct solar irradiances and of atmospheric emittance, Data acquisition 
was computer automated in February 1982 and, since then, one-minute inr>egrated 
radiometric data, as well as one-hour integrated data, have been recorded. 

Work has continued on the correlation of global insolation with cloud 
cover fractions using the ARL model, The results for the first complete 
year of data are discussed in part Ill-C of this report. A theoretical model 
has been developed and a data set has been accumulated to test this model 
which is used to estimate global insolation from satellite radiometric data, 
Detailed information about the model is provided in section IV of this report. 
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I. SOLAR RADIATION MEASUREMENT PROGRAM 

An observation platform for solar and atmospheric radiation measure- 
ments was established on the roof of Turner Hall (latitude 37.02°N, longitude 

¥ 

76,34^, and elevation 24 meters) February IS, 1981, A radio tower and a 
smokestack are the only two obstructions greater than ten degree^l above the 
horizon. Routine measurements are made of global solar irradi^ince, diffuse 
solar irradiance, direct solar irrsdlance, and atmospheric emittance. Infor- 
mation (Griffin, 1982) about the solar radiation measurement program was / 
presented at the April 20-23, 1982 meeting of the Virginia Academy of Science 
held in Blacksburg, VA. The Abstract of the presentation by T, J. Griffin is 
attached as Appendix I, A summary of the solar and atmospherlit data available 
at the time of this report Is provided In TABLE I. 

A. Ins trumenta tlon 

Since detailed descriptions of the radiometric Instrumentation are avail- 
able In the first Annual Report (Whitney, 1981) for this grant, only a summary 
of our measurement capabilities Is presented here. The global solar irradiance 
on a horizontal surface is measured by an Eppley Precision Spectral Pyranometer 
(PSP) with a hemispherical WG7 Schott glass dome. Diffuse solar irradiance on 
a horizontal surface Is measured by an Eppley PSP with a WG7 dome that Is mount- 
ed on an Eppley Solar Tracker and Occulting Disk System. Direct solar irradiance 
is measured by an Eppley Normal Incidence Pyrhelioroeter (NIP) mounted on e Solar 
Tracker. An Eppley Hickey-Frieden Absolute Cavity Pyrheliometer Is used to 
calibrate the NIP regularly. Atmospheric emittance Is measured with an Eppley 
Precision Infrared Radiometer (PIR). In addition, turbidity (Volz, 1974) 
measurements are made approximately hourly for clear sky conditions using a 
VOlz Sunpho tome ter. 
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The all-eky camera system was modified in order to allow computer 
activation of the camera. The Aetna fiaheye adaptor lena and previoua Ricoh 
XR-2 lena syatem were mounted on a Fentex ME- II camera body with an electronic 
shutter and an autowind ayatem. The camera is mounted on lin Eppley Shadow 
Band Stand and aimed a t the tenith in order to eliminate the need for frequent 
adjustment of the sunshade. A clock and a date card are placed on the inaide 
face of the ahadow band to document the time of day and date that each photo- 
graph is taken. Whole-sky photographs (color slides) are taken every one-half 
hour on weekdays when clouds are present without precipitation. Color slides 
are uaed in order to aid in the diatinction between clear sky patches and dark 
cloud bottoms. The image of the sky is projected directly onto an analysis 
grid to obtain the fractional cloud cover. Processing and analysis tlmea for 
whole- sky photographs have been greatly reduced from the previous black-and- 
white enlargement method. 

The wavelength range of each of these instruments is listed in TABLE II. 
Information about the measurement frequency and about the time periods that 
Insolation data are available is presented in TABLE I. 

B . Calibration 

Bach radiometer was calibrated prior to delivery by the manufacturer. 
Subsequent calibrations have been performed by comparison to secontjiary or 
primary standards at Hampton Institute or at Eppley Laboratory. The NIP was 
calibrated by comparison to the Hickey-prieden Absolute Cavity Pyrheliometer 
which is considered to be a primary standard. The pyrgeometOr (PIR) and the 
older pyranometer (PSP # 20022F3) were recently recalibrated at Eppley Lab- 
oratory. The other pyranometer was then recalibrated by comparison with the 
calibrated pyranometer by operating side-by-side for a period of two days and 
comparing the two-day integrated global irradiance totals. Calibration data 
are presented in TABLE III by listing each Instrument , the date and site of 
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each calibration, the aeneitivity factor, and the percent change in each 

initrumenUVa eeniitivity factor, Comparieono between the two pyraiibaietere 

"ii 0 

indicate a conaiatent difference of about one percent which ie within the 
two percent accurecy of each inetrument. The rather large change in oeneitivity 
of the older inetrument implies either that both inetruments have changed 
sensitivity at the same rate, or that an error exiets in the 8 - 301^62 calibra- 
tion. Discussions with Bppley Laboratory about this matter ere currently be- 
ing pursued and correction factors for each instrument will be calculated for 
each month of operation since the sensitivity factors for these two pyrenometers 
apparently chsnge,d by more than two percent. A change In calibration standards 
St Eppley Laboratory in October 1961 accounts for 2.5% of the sensitivity factor 
changes. The second pyrsnometer is being sent beck to Eppley Laboratory for 
recalibration to verify this apparent change in sensitivity factors. 

The recorder systems were calibrated at least every six months by using 
a stable millivolt source and adjusting the integrator and strip chart recorder 
gains to give the proper readings, The maximum adjustment made for any of the 
integrators has been 0.4% while larger adjustments were occasionally, needed 
for the recorders. The radiometric data are obtained from the integrator LED 
display BCD panels and are only rarely supplemented from the strip chart 
recorder trace. 

The Volz Sunphotometer is calibrated by checking for linearity of the 
sensor on a Langley Plot of readings versus air mass, 

G. Meteorological Data 

Standard hourly meteorological observations taken at nearby Langley Air 
Force Base (LAFB) are picked up on a monthly basis from Detachement 7, Third 
Weather Squadron. These data include information about cloud height, cloud 
type, fractional sky cover, precipitation, sea level pressure, and temperature. 

These data are supplemented by whole-sky photographs and visual cloud obser- 


votlons at Hampton Institute. Additional ueteorologioal data and ozone data 
are purchased from the National Climatic Center archives in Asheville, North 
Carolina and The World Ozone Center, Environment Canada, Toronto, Canada for 
use in the data analysis. 

D . Data Handling and Quality Control 

A data storage procedure £or the radiation data has been devised to effic- 
iently handle the data and ensure quality control. First, the integrated radio- 
metric data and times are stored on a Tektronix 4051 microcomputer* a internal 
magnetic tape unit. Then on a weekly basis these raw data files ere transferred 
via computer hookup from the Tektronix 4051 to the PDF 11/34 minicomputer where 
the data are permanently stored on 1600 bpi magnetic tape. All preliminary 
processing is done on the Tektronix 4051, while data analysis routines and 
application programs are performed on the PDF 11/34 system. The data are 
examined for errors by a computer program that locates the gaps in the data and 
Identifies errors to be corrected. 

The automated data acquisition system for global and direct solar irradi- 
ances was fully operational February 1, 1982, and since then, integrated radio- 
metric data for one-minute intervals have been obtained directly from the 
Bppley Integrators and stored on magnetic tape. Specific data handling proce- 
dures and quality control for these data are discussed in more detail in the next 
section of this report. 

Prior CO the installation of the automated data acquisition system, only 
one-hour Integrated measurements were recorded on magnetic tape. Printed data 
from the integrators were scanned on s daily basis for missing or problem data 
and incorrect timing caused by power failures or other electrical and mechanical 
malfunctions. Hissing data were supplemented by the strip Chart record when 
available. On a weekly basis, these data were manually entered into ’the 
Tektronix 4051 microcomputer, inspected for operator errors, and then transferred 
to the PDP 11/34 data storage tape. 
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II. AUIOHATBD DATA ACQUISITIOH SYSTEM 

One of Che major objectlvea of Che second year of Chis research projacc 
was Co develop and ImplemenC an aucomaced daCa acquis ICion sysCem. This sys- 
Cem reduces Che manpower requiremencs . and Che number of operaCor induced 
errors, involved wich daCa cransferral tt<m Che InCegra Cor- recording sysCem 
Co final permanenc magneCic Cape sCoraga. The hardware requiremenCs , operacing 
procedures, feaCures, and performance of chc sysCem are discussed below. 

A. Svseem Overview 

The auComaCed daCa acquisiCion sysCem conneeCs Che radiomeCtic sensor- 
inCegraCor recording sysCem wiCh Che rolcrocompuCer I/O capabilicies via an 
inCerface box. The microcompuCer reads Che inCegraCed radioroeCric values 
from Che BCD inCerface panels ac a presec Cime inCerval and records chese 
daCa on Che inCernal magneCic cape unic. Whole~aky phoCographs are criggered 
by computer command aC preset times and a record of Cime and picture number 
is recorded on Che magneCic tape. Integrated atmospheric emiCCance, and 
global, direct, and diffuse solar irradiances are measured at one-minute 
intervals from 0400 EST CO 2000 EST and at Cen-minuce intervals for Che rest 
of the night. The data for each month are transferred Co separate 1600 bpi 
magneCic tape for permanenc storage. 

B. SvsCem PescriDtion 

The data acquisition system can be separated into three subsystems: 

(1) the radiometric sensor- integra tor system; 

(2) Che integrated data sampling and recording system; and 

(3) the data storage sysCem. 

A flow chart of the data acquisition system is presented in Figure 1. The 
radiometric sensor-integrator system is located at the top half of the chart 
and includes the following five components: 

(1) the Eppley radiometers; 
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(2) th« BppUy intagratori that aoa tha inatantanaoiia raadinga le<M tha 
radi,oaMtara; 

(3) tha Digitac printara th*t provida a hard copy printout oi ona-hour 
radlanca valuaa; 

(A) tha X-Y atrip chart racordar.s]that provida a hard copy of inatan* 
tanaoua radiancaa; andi 

(5) tha camara ayataia that takaa tha vhola-aky photographa. 

Tha intagratad data aampling and racording ayataa locatad at tha lowar laft-hand 
aida of tha chart includaa: 

(1) tha intagratoi' aignala that ara providad at tha BCD intarfaca on tha 
Bpplay intagratori; 

<2) tha intarfaca box that cantralizaa tha data for computar accaaa; 

(3) tha ROMb that intarfaca tha BCD data to tha Tektronix 4051; 

(4) tha Tektronix 405L microcomputar that raada the ROMa and aetivataa tha 
camara photoeounter; and, 

(5) the magnetic tape cartridge where tha computar atorea tha data. 

The data atoraga ayatem ia located at the lower ri^ht-hand aida of the chart 
and conaiata of three componenta: 

(1) another Tektronix 4051 microcomputer with a RS232 intarfaca that acta 
aa a link between the PDF minicomputer and the internal tape of the 
Tektronix computer; 

(2) the PDF 11/34 minicomputer ayatem that reada the data from the 
Tektronix 4051; and 

(3) the 1600 bpi, 9 track magnetic tape for permanent atorage of the data. 
The hardware uaed in the integrated data acquiaition and atorage ayatema are 
Uated in TABLE IV along with two other devicea that were uaed In the development^, 
and teating of the Interface box and aoftware programa. Detailed information 
about the operation of these devicea and about the construction of the interface 
box la kept in a documentation file in the Solar Energy Measurement Laboratory. 

C . Operating Procedure 

The data acquiaition system operates on the Tektronix 4051 microcomputer 
using the computer program SOLAUTO written in BASIC by D, D. Venable, The 
program works by comparison of the time provided by the Real Time Clock ROM 


7 


»» 




Pack wlfch Cimaii that ate calculated from fixed time Intervila entered into the 
program. At the proper tirnea the computer reada each ROM for the radiometric 
data, or activatea the camera ahutter and autowind ayatem. A flow chart for 
the SOLAUTO program la praaented in Figure 2 in order to repreaent the logic 
uaed in thia part of the program. A copy of the program liating ia kept in 
the documentation file in the Solar Energy Maaaurement liaboratory. 

The correct time, date, aelected aampling tlmea for the radiometric 
data and for camera activation, and alao, uaer meaaagea can be entered into 
the computer by uaing the Uaer Defined Keya (UDK) on the computer. The radio* 
metric data are recorded on the magnetic tape in 72 character recorda. The 
firat two-digit flag (0#) of thl> data record indicatea the type of data atring 
uaing: 01 for radiometric data, 02 for a photograph record, 03 for a uaer 

meaaage, and 04 for a ayatem meaaage. The next charactera indicate the day of 
the week (three letter abbreviation) aeparated by two apacea from the date in 
the form: day-month-year (DD-MMM-YY) . The date ia repreaented by two diglta 

each for the day and for the year, and by the three letter abbreviation for the 
month. The next aet of numbera (HH:MM:SS) in the data atring are the time 
(EST) repreaented by two diglta each for the hour, minute, and aecond, and the 
remaining charactera are either a meaaage or data. A radiometric data record 
ha a the form: 

01 , DDD DD-MMM- YY HH: ^fl: SS , GL0 , NNNNN, DF# , NNNNN, IK# , NNNNN , DR# , NNMHN 
The radiometric data are recorded aa integrated totala atarting from aero 
at midnight uaing the two letter abbreviation for each radiation component, 
a single digit activation code (1 for on, 0 for off), and five digits for each 
reading (NNNNN). Global solar irradiance is abbreviated by GL: diffuse solar 
by DF: infrared (atmospheric emittance) by IR: and direct solar by DR. 

The system is Initially operated by loading the SOtdVUTO program into the 
memory of the Tektronix 4051 computer and then inserting the data tape cartridge. 
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A£tftr enttring "RUN"t th* p«ranMit«ri art §%t to indleito whi^h RONa ahould b« 
read» tht fraquoncy for aampllng data and activating tha caMrii tha propar data 
and tina, and tha data atoraga/diaplay looationa* TABLE V Liata tha program 
faaturaa by Idantifying tha function of acch kay, whita tha Hoat fraquantly uaad 
kaya ara diacuaaad balov. UDK 1 - *'Nota to Tapo'' ia uaad to aiaka coaaaanta 
•bout lenaor cLaaning, tape changaaj, or othar inportant naaaagaa, UDK 2 ’’Raad 
ROMa" li uaad to chaok tha R(M raadinga by coapariaon with tha LED diaploya on 
tha integratora, Tha caowra ia manually activatad by UDK 3 ~ '^Photograph'.' and 
la fraquantly uaad to chack tha alactronic ahuttar oparation and to varlfy tha 
photograph numbar. The photograph or aanaor reading time Intarvala ara changed 
by ualng UDK 5 <• "Intarvala" and tha Photograph numbar can be raaat or corrected 
by ualng UDK 10 - "Photocount". In order to Interrupt the operation of any 
Inatrument for maintenance or calibration without cauaing a fatal program arror^ 
the correapondlng ROM can be turned off (or on) by uaing UDK 7. Program 
execution can be atopped completely by ualng UDK 20 <* "Stop" and raaumad by 
using UDK 19 - "Continue". The use of UDK 13 "Status" provides a quick check 
of the I/O status, aampllng rates, photograph numbers, the date, and the time 
that are currently stored In memory. 

In addition, several operations can be performed with the data during the 
sampling operation. Because of the time required to rewind and read the mag- 
netic tape cartridge, the data atrlnga are internally stored by the 4051 computer 
when using UDK 15, 16, or 17 In order hot to lose any data. The function of UDK 
15 - "Tape Dump" ia to read the tape and display the data for a certain day, or to 
list all photograph times, user messages, or system meassges. By using UDK 17 - 
"Graph" a plot of global, direct. Infrared, or diffuse Irradlancea can be obtained 
on the Tektronix 4051 Graphics Display Screen. New data tape cartridges are 
installed approximately once a week. 
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Th« p«rformanc« o£ fcht eonpuuif tutooiatiad data acqulaltiQii cm be 
nmaured by calculation oi thi anMunt ot data loat in coiapatiaon with tha data 
racovarad on Mgnatic tape* Tha parformanca racord o£ tha pcaviduo data 
acquiaition ayatam muat alao ba conaidarad ainca both ayatama have aoaa of 
aaaa cauaaa for loaa of data. Pot axamplaf diffuaa colar radiooiatrid data ara 
not collactad during calibration parioda for aithar pyranonatar. Intarcooipariaon 
of tha two pyranomatara ara mada for tha horizontal global aolar oriantationi 
which maana that no diffuaa aolar data ara availabla during tha calibration 
parioda. Elactrical atorma can inter fara with computer program axacution and 
can atop data collection until the ayatem ia raatored to normal operation. 
Electrical atorma, which alao affect the printera and intagratora, reaat the 
printer timea and the integrator count valuea, an effect that deatroya the 
printer data during the time period of the atorro. The atrip chart recordara 
can be uaad to retrieve most of the loat data and have been uaed to obtain 
hourly Integrated data. Mechanical solar tracking failures can destroy the 
direct or diffuse aolar data obtained by either system. 

One*’hourly integrated data availability ia presented by month in the DATA 
RECOVERY RECORD TABLES (VI- IX) for the full period of the insolation measurement 
program at Hampton Institute. TABLE X contains the one and ten minute integrated 
data recovery information specifically related to the automated data acquiaition 
syatem. The measurement interval for global and direct solar irradiances was 
set at one-minute starting at 0942 EST February 1, 1982. One-minute integrated 
data sampling for the diffuse Solar and atmospheric emittance began at 2038 EST, 
March 1, 1982 when the second half of the interface box was completed. On March 
13, 1982 the sampling Int#' ^al at night was changed to ten minutes (from 2000 
EST to 0400 EST) in order to reduce computer tape storage requirements. One 
Tektronix Data Tape Cartridge is used to store approximately six and one-half 
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(6^) dayi of Inaolation data. 

maaauraa of tha performanca of tha automatad data acquialtion ayatam 
ara to calcuLata tha ayaraga tlma parlod batwaan falluraa and tha avaraga 
langth of tlma loat for aach error. Thaaa maaauraa are rafarrad to a# Mean- 
time to Pall and Downtime per Error, raapectlvely. The Mean Time to Fa)Vi can 
be calculated by dividing the total poaalbla number of data racorda by tha 
numbar of falluraa and by the data record aampllng rata, aa long aa the aamp« 
ling rata la conatant. The Downtime per Error can be calculated by rdlviding 
the number of data records lost by the number of errors and by tha aampllng 
rate. The results of these measures of performance ware reported (Blakay, 
1982) at the 39th Joint Annual Meeting of the National Institute of Science 
and Beta Kappa Chi Scientific Society held In Washington, DC, March 17'20, 
1982. Rody Blakey, an undergraduate assistant on this project, used the data 
obtained during the first month of computer automated data acquialtion to 
calculate the Meantime to Fail and Downtime per Error. The results of his 
analysis were an everage of 8.7 days between failures and an average of 6.7 
hours downtime per error. The large downtime per error waa caused by the two 
nlghtlme failures that stopped program execution for several hours until the 
system could be returned to normal the following working day. (Mr. Blakey 
received a Third Place Award for his paper in the Mathematics and Computer 
Science Section of the Meeting and, also, one of two general awards given.) 

A copv of the abstract of hla paper la attached to this report as Appendix II. 


11 


M 




III. RADIOMETRIC AND METEOROLOGICAL DATA ANALYSIS 

Several typea o£ data analyals have been completed using the radio- 
metric and metaorological data collected ac Hampton Institute., The first 
type of analysis is to present average values for the various measurements 
in graphical form. A second type of analysis is required to treat the raw 
turbidity measurements in order to obtain useful parameters such as optical 
depths and preclpitable water. Correlation of data with empirical formulae 
is a third type of analysis performed. At this time, satellite-derived 
radiometric data have not been analyzed by any of these methods but are 
expected to be completed within the next few months. The results of the 
analysis methods used are discussed below. 

A. Radiometric Data Presentation 

Diurnal variability in the various solar insolation components are obser- 
ved by comparison of the average hourly values for afternoons with values for 
mornings. Average hourly global, diffuse, and direct solar irradiances and 
a tmospheric emittance are plotted for each month from July 1981 through 
June 1982 in Figure 3 through Figure 14. Diurnal variability is indicated 
by the lack of symmetry in these graphs. Seasonal variability can be seen 
by plotting global solar irradiance for clear sky days selected from each 
season such as is done in Figure 15. Plots can also be made for values 
averaged over shorter time intervals (one minute to sixty minutes). Five 
minute integrated values are plotted routinely for each day to verify the 
computer system data. 

B. Turbidity Data Results 

The Volz Sunphotometer is used to obtain atmospheric turbidity para- 
meters a and $ , the aerosol extinction coefficients or optical depths at 500, 

875, and 390 nanometers, and preclpitable water concentrations. The Instru- 
ment operates by generating a dc voltage proportional to the direct normal 
solar Irradiance measured through narrow band transmission filters centered 
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on 390» SOO, 875, and 945 nm, Measucements ace made under cloud free sky * 

conditions. Calculations are then required to convert voltages into the 
turbidity parameters. 

Aerosol optical depths at the wavelengths measured are calculated using 
the following equation from Volz (1974): \ 

^ In ( F Iq/I ) - K ( + Tq ) ; where, 

aerosol optical depth, 

Iq" extraterrestrial irradiance at mean polar distance, 

I ■ irradiance measured with the Sunpho tome ter, 

F p correction factor for mean sun-earth distance 
H optical air mass, 

Rayleigh optical depth, 

Tq" ozone optical depth, and 

K ■ correction for differing optical path at M > 2.5. 

Monthly average aerosol extinction coefficients at 390, 500, and 875 nm 
for the measurement period March 1981 through June 1982 are graphed in Figures 
16-18, respectively. Turbidity parameters 0 and 8 are calculated from the 
aerosol extinction coefficients using the Angstrom (1961) empirical expression: 

Ta Cll ■ ; where, 

1 «• the wavelength in mlcrone, 

Ta ClJ. “ the aerosol optical depth at wavelength 1, 

0 ■ RngstrSm turbidity coefficient, and 
a » XngstrWra wavelength coefficient. 

The turbidity coefficient 8 is a measure of the quantity of aerosol 
particles in the atmosphere. Large valuC'S of 8 indicates that a high aerosol 
concentration exists in the local air ttass. Small 8 values indicate a clear ' 

air mass with low aerosol concentrations. For the measurement period, the 
largest clear sky daily average 8 value was 0.48 obtained on July 28, 1981 
when atmospheric conditions were extremely hazy. The smallest daily average 
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was 0.05, obtained on February 28, 1981. 

Fl8u>^< 19 graphs monthly average B values for the measurement period. 

The graph shows that local aerosol concentrations of dust, smoke, sea salts, 
and other particles peak In the midsummer. 

The Xngstrbm exponent a Is a measure of aerosol slse distribution. 

A large a value Indicates small particles dominate the aerosol populatloni 
while a small a means that large particles dominate. Xngstrbm states that a 
typical value for c* Is 1,3. At Hampton Institute we have obtained a range 
of dally average « values from 1,73 to - .89, with an average Value of 0,72, 
Figure 20 plots monthly average values of ot and, as shown by the graph, 
the values are highly variable. This variability may be due to the fact that 
local aerosol distribution Is a function of the type of air mass over the 
measurement site. Further study Into the variations of a will focus on the 
particle distribution under different air mass conditions. 

C, Correlations with Cloud Cover Fractions 

An empirical model developed by the Air Resources Laboratory (ARL) of the 
National Oceanic and Atmospheric Administration was Initially selected for 
correlation of the global solar irradlance data with fractional cloud cover 
values. This Is an established model used by meteorological stations, and 
Is used to correlate: hourly average global solar Irradlance with opaque cloud 
cover fractions which makes It particularly compatible with our solar energy 
and cloud cover measurement program. The results of the correlation can be 
used for comparison with other measurement sites In order to characterize 
the measurement site at Hampton Institute. 

The major factors that affect the amount of global solar energy received 
at the ground on a flat surface are solar zenith angle and opaque cloud cover. 
The ARL empirical equations selected for use with our measurements are: 

(1) SRC - Aq + Aj^ cos ZA + Ag cos ^ ZA + Aj coS ZA, and 
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(2) SR - SRC <Bq + OPQ + B 2 OPQ^ + B 3 OPQ^ + B^ RN) 

SRC If thf folft rfldlation hourly value for clear aky oondltiona; 

ZA ia the aenlth angle at the midpoint of each one-hour interval: 

OPQ ia the average opaque cloud cover fraction; and 

!f<N ia a rain term that ia equal to one if aome form of precipitation 
la reported, othetwiae it la zero* 

The coefflcienta for the firat regreaaion equation are calculated 
aeparately for mornings, and afternoons, for each month of the year in order 
to partially account for diurnal and seasonal changes in atmospheric tq^Thidity, 
water vapor, and other such factors. The first and last partial hours of the 
day are not included in the ttcgreasion calculation. The coefficients for the 
second equation tire normally calculated for mornings plus afternoons using a 

j 

full year of data at a time, 

1. Clear Sky Model Results 

The clear sky data have to be analyzed before the cloudy sky data can be 
normalized by the solar zenith angle dependence. Determination of clear sky 
hours was made using the LAFB cloud observation data set. A plot of hourly 
values of clear sky global solar for mornings and for afternoons are shown in 
Figures 21 to 33 for the months of June 1981 through June 1982. The line 
plotted through the data was obtained by the Marquardt's nonlinear least- 
squares method (Marquardt, 1963) using Equation 1. The regression coefficients 
and number of points used In each correlation are listed for each month in 
table XI. The number of tot«Ally clear sky hours was insufficient for meaningful 
analysis for most months and, therefore, were supplemented with ''nearly clear" 
points for which the cloud cover fraction was less than 0.2 and the strip chart 
trace showed no indication of clouds. These nearly clear points are added in 
order to increase the number of points in the data set used in the regression, and 
to extend the meaningful range of the curve fit to zenith angles for which 
totally clear sky points are not available. Nearly clear sky data points were 
selected from hours which were coded as clear at the beginning or end of the 
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hour, and alao from houra which war# coded aa having low frachlona of 
tranaparent or aamltranaparant clouda^ In Figure 21 ho 33 the clear 

■'ll 

aky data polnta are plotted with open clrclea and the neariyt^claar aky 

!i 

H 

polnta with plua algna. Reaulta for the data prior to June 1981 were, 
reported In the flrat Annual Report (Whitney, 1981) 

2, Cloudy Sky ARL Model Reaulta' 

The regreaalon coefflclenta for clear aky mornlnga and for clear aky 
afternoona.were uaed In Equation 1 to eatlmate the clear aky Irradlance for 
each one>hour Interval. The cloudy sky data were then normalised and fit 
to Equation 2 by using the Marquardt's nonlinear least-squares method. The 
dftta for mornings and afternoons are plotted for each month from June 1981 
to June 1982 In figures 34 to 46. The cloud cover fractions obtained from 
whole sky photographs are represented in these plots by plus signs and 
satellite photoprint-derived fractions are ropresented by open circles. 

The correlation coefficients for each month of combined morning and 
afternoon data are listed in TABLE XIl. The cloud cover fractions used in 
these correlations are averages of the LAFB visual observations made at the 
beginning and end of each hour. Coefficients for data sets with more than 
230 points, the limiting number of data points in the analysis program, were 
obtained by calculating the average normalized irradlance for clear sky 
conditions and for overcast sky with and without precipitation. The first 
coefficient Bq was assigned the mean SR/SRC value for clear and nearly clear 
sky conditions since Bq Is the value of Equation 2 at the point 0PQ“0, i. e. 
for clear sky. The fifth coefficient B^ was assigned the mean difference 
between values with precipitation and values without precipitation at total 
cloud cover, since the rain term RH nearly always equals one for overcast sky 
conditions (0PQ»1) and is rarely non-zero for partly-cloudy sky conditions. 

ARL coefficients for Equation 2 are normally reported for a complete 
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calendar year of data. The reaulta for the first year of data,. Match 1981 
through February 1982 are presented in tABLE XIII for both the olear-sky 
and cloudy-sky equations. Since the cloud fractions obtained from the LAFB 
visual observations were too numerous for the computer program usedi the 
average normalized global irradiance values (SR/SRC) were calculated for each 
0.05 cloud fraction step. The mean values and points one standard deviation 
above and below these means were used to calculate the correlation coefficients. 
All the data points for partly-cloudy sky with precipitation were also included 
in the calculation. The correlations were repeated with fractions obtained 
from satellite photoprints and with the whole sky photograph values using 
the mean values for clear, overcast without precipitation, and overcast with 
precipitation for the endpoints. The results are given in TABLE XIII and 
the data are plotted in Figures 47 and 48. 

The large scatter of the data in Figures 47 and 48 imply that the cor- 
relation coefficients are not defined well enough to draw many meaningful 
conclusions. The results for the cloudy sky model ace similar to those from 
Other meteorological recording stations along the Bast Coast but the clear 
sky coefficients are not. The probable error in each coefficient is calculated 
in the analysis program used and the errors in A^ and A^ for the clear sky data 
are frequently much larger that the coefficients themselves. Thus a comparison 
between the measurement site at Hampton Institute with other sites wf^uld 
require a much larger data set before the ARL model could distinguish between 
sites with similar climates. The scatter in the data indicates that fractional 
cloud cover and precipitation are not the only parameters chat affect hourly 
global insolation at a given location, Additional factors involving cloud 
geometry, locations, altitude, and thickness are significant. Perhaps, the 
method of calculating fractional cloud cover is flawed, and a single fraction 
for an hour should be calculated from some type of time averaged set of 
fractions. These fractions should not be just the cloud cover fraction directly 
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over the eite, but retheti the fteotion between the ilte and the aun« 
Normallziid global irradiaiiGe ia plotted veraua whole-aky photo fractlona 
in Figure 49 and veraua Satellite photoprint fraotiona in Figure 50* The 
error bars indicate the mean values for overcast sky, plus/minua one atandard 
deviation. 

The frequency of occurrence for certain average cloud cover fractions ia 
plotted in Figures 51, 52, and 53, from the three different cloud cover 
fraction sources used in the correlations with the ARL model. There were 
2,926 fractions used from the LAFB visual observations, 53 fractions from the 
satellite photoprint analysts, and 95 fractions from the whole-shy photograph 
analysis. The 0.05 interval in cloud fraction is the result of averaging «two 
fractions provided in tenths from the LAFB data set, and from rounding the 
fractions from the other data sets to the nearest 0,05. The actual number of 
whole-sky photographs is much larger than 95 since only hourly averages were 
used. 

The normal blmodal distribution of cloud cover fractions is indicated in 
Figure 49, but is not observed in Figures 50 and 51 since photoprints and 
photographs were not used for either clear sky or overcast sky conditions. 

The slight peak at small fractions for the satellite photoprint data set 
is probably due to the inability to distinguish very low If/Vel and very small 
clouds from the ground in the photoprints. The whole-sky photograph data eet 
indicatea that the cloud cover fractions for partly'-cloudy sky are distributed 
uniformly for the small region of view Cof the order of a few kilometers 
in radius) as well as for the much larger region used In the LAFB visual 
observations (horizon to horizon). 
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IV. theobetioal modeling 

the theoretical model developed by J* R, Foremen, waa aeleoted for 
uae with our data and thia method is preaented be'loiir in detail. The clear 
aky part of thia model ia being applied to a March 1979, Michigan regional 
data aet aa a PhD, in Atmoapheric Science diaaertation project. Our local 
data aet ia being uaed to teat and develop the model on a amall acale and 
to expand the model to partly cloudy aky conditions. The month of February 
1982 was selected as the flrat month to apply this model for the following 
reasons: (1) the ground albedo was expected to remain nearly constant 

throughout the month providing that there was no snow and that there was 
little change in the amount of ice in the tidal basin (2) the automated 
data acquisition system provided short-time interval (one minute) integrated 
data; and (3) clear sky conditions were sufficiently common to provide a 
data base for characterization of ground albedo and normal atmospheric absorp- 
tion and scattering. 

A • Summary of the Model 

The method of estimation of global solar irradiance on a horizontal 
surface is an extension of the short wavelength energy balance equation 
developed by Ellis and Vender Haar (1978) 

^bg - jdo - Ir " Ifl) where; 

^hg equals the horizontal global shortwave irradiance at the earth's 
surface; 

a equals the local ground albedo; 

^o equals the horizontal extraterrestrial solar irradiance weighted 
by the spectral response of the satellite imaging device; 

^r equals the total shortwave irradiance reflected to space by the 
earth's atmospheric system; 

and, ^a equals the total (incident and reflected) shortwave radiation 
absorbed by the various components of the atmosphere. 
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Thfr tonq ii oalcuUtcd £rom the ■atellite meeeured brlghtneee of a 
line element and an appropriate angular reflectance model, which la being 
■elected from the recent lurvey by Davie and Cox (1981), 

The portions of due to absorption by ozone and water vapor can be 
found to high accuracy from parameterizationa developed by Laoie and Hanzen 
(1974) and the absorption by permanent gazaea can be calculated using a model 
presented by Yamamoto (1962). The absorption of light by dry and semi-wet 

H 

aerosols will be parameterized to the mean low level relative humidity, f, 
as found from a weighted average of the mean surface to 500 mb, the 850 mb, 
and the 700 mb relative humidities. The weighting will be done in accordance 
with some mean vertical distribution of aerosols in the low levels, according 
to a formula suggested by Hifnel (1976): ■ I^^exp{ -Cjm(^) (1 - f)^2 ); 

where is the total radiation absorbed by atmospheric aerosols; 

loj is the incident radiation with ozone absorption accounted for; 
p is the unreduced surface pressure (p^ equals 1013.25 mb); 
m is the relative optical air mass of Kasten (1966); and 
C| and C 2 are empirical coefficients. 

The fourth absorption portion of is caused by clouds and is the least 
understood of the various absorption mechanisms. Several relationships between 
absorption by clouds and the total short-wave radiation reflected upward by 
the clouds, the presence or absence of precipitation in the clouds (as revealed 
by radar reports) and the presence or absence of ice crystals in the cloud 
top (as inferred from the infrared estimated cloud top temperature) will be 
investigated during the next year. 

B. Data Sets 

The meteorological data needed for this model have been purchased from 
various sources. Most of the airways reports and weather data were purchased 
through the Environmental Data Services, National Climatic Center, Asheville, 
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NC, Some o£ these data are avalleble In August and were not received at 
the time o£ the report. The local study region is indicated in Figure 54. 
by the dotted lines and with the Hampton Institute measurement cite indi- 
cated by an X. As noted in the legend, each meteorological site is repre- 
sented by a separate number. The sources o£ data outside the study region 
are used to establish parameters at the boundaries o£ the study region and 
as back-up data sets since several data sources are not regular one-hourly 
reporting stations. A summary o£ the data sets needed in the analysis is 
given below in outline £orm. 


1. Satellite Data 


a. Source: GOES-East visible and in£rared digitized brightness values. 


b. Resolution; 0.9 km x 0.9 km visible, 3.7 km x 3.7 km in£rared. 


c. Data Array Coordinates: 

Point 
Center 
NW Corner * 
HE Corner - 
SE Corner - 
SW Corner - 


Satellite 
line X element 
2992 X 7517 
2942 X 7567 
3042 X 7567 
3042 X 7467 
2942 X 7467 


Geoaraphical 
latitude x longitude 
37.019°N X 76.338°W 
37.587°N X 76.884°W 
37.584°N X 75.80l°W 
36.470°N X 75.798°W 
36.473°N X 76.858°W 


d. Period; Clear Sky Hours for Feb. 7,8,l0jll,22 and 28, 1982; 

Cloudy Sky Hours for Feb. 4,13,17,18,21, and 24, 1982; 

2. Meteorological 


a. Insolation; 

I. Type; hprizontal global - 0.3 pm to 2.8 pm. 

II. Frequency; ten-minute totals centered on satellite scan time (GMT) . 
ill. Source; Eppley PSP one-minute integrated data at Hampton Institute. 

b. Humidity 

i. Hourly surface airways reports of temperature, dew point, and 
altimeter settings from local airport, Coast Guard, Air Force, 

Navy, and Army weather data sources for low level data. 

ii. Radiosonde data from Wallops Island and Sterling VA and from 
Cape Hatteras and Greensboro NC for upper level data, 

c. Cloud Parameters 

I. Cloud fraction: HI photographs and satellite brightness readings. 

II. Cloud type and height; Langley Air Force Base observations. 

HI. Water content: radar reports (facsimile charts) from local 

airports. 

iv« Cloud top temperature; Infrared satellite data. 
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d« Ozone end Aiitoiole 

i* Tuvbldity fiidingi eb Hempton Inibitute 
li, Dobeon epee fcropho tome be r reading! of total ozone column over 
Naahville, TH/ 'rallahaaiee, FL, Wallop! l!land, VA, ond 
Wa!hington> DC. obtained from the World Ozone Data Center, 
Environment Canada, Toronto, Canada. 

3 . Miicellaneou! 

a. Ground albedo; calculated from clear fliy iatelllte brightneae 
value! ueing an equation developed by Otterman and othera (Otterman, 
1980) . 

b. Elevation; moat points within the study region are treated aa being 
at aea level but uae of United States Geological Survey Maps is 
being made to ascertain elevations throughout the study region. 

c. Cloud shadows; use of photoprint and calculation of shadow positions 
from sun-cloud-satellite geometery, 

C, Progress 

The first group of steps in the theoretical modeling have been completed; 
i.e., a model has been selected and a data set has been accumulated for com- 
parison with the model. Satellite radiometric Imagery has been obtained from 
the Environmental Data Services of HOAA for the Hampton Institute study region 
for selected times within the month of February 1982. Meteorological data 
have been ordered from the National Climatic Center for local observation 
sites for the same time period. The ozone data has already been obtained 
and the horizontal global irradlance values are readily available. The pro- 
cedure for the application of J. R. Foreman's model is outlined below to 
give an indication of the steps that must be performed* 

1. Assemble Data Set 

a) select times and dates; 

b) order satellite imagery (visible and infrared brightness values); 

c) request meteorological data (ozone and relative humidity data); 

d) calculate satellite scan times for insolation intervals; 

e) assemble horizontal global insolation values; 

f) calculate relative humidity for surface to 500 mb, 700 mb, and 850 mb 

g) divide data set into clear and cloudy data sets. 


m 
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2. Uft Qltar Difc« Stt tot 

•) caloulitt ground olbido (ct) from IL » I nnd anguUr roflootanot modtli 

b) . ch«r«cUrirt •baorpfclon by pttmintnfc §•■•••;[ 

c) colculato «btorpi:ion by oson* from ipiorrophotoRMfctr or VoU 
■unphofcomtfctr data; 

d) obtain coafficlanta for abaorptlon by aaroaola from 1^, {» and 
abaorption by oaona and parmanant gaaaaa, 

3, uaa Cloudy Data Sat to» 

a) aoparata cloud ef facta from clear aky af facts; 

b) datermlna dipandonci on cloud typS) halghti transmlttanca and amount. 

Since the data sat has not been completely assamblad» we have not pro- 
greiaed to the ground albedo and model parametar calcttlitlonf » Ihe compleca 
data set should be assembled within the next month and then work can commenee 
on the application of the model. 


» t I 
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TABLE II 

RADIOMETRIC INSTRUMENTATION AND WAVELENGTH CAPABILITIES 


Inatrument 

Eppley Precision Spectral pyranometer 

Eppley Precision Infrared Radiometer 

Eppley Normal Incidence Pyrheliometer 

Eppley Hlckey-Frleden Absolute Cavity 
Pyrheliometer 

Additional Wavelength Ranges for the 
Pyranometer s and Pyrheliometer s 

Volz Sunphotometer 


Normal Observation Wavelength Range 

0.285 to 2.8 microns 
4.0 to 50.0 microns 
0.285 to 4.5 microns 
0.2 to 50 microns 

0.53 to 2.8 microns 
0.63 to 2.8 microns 
0.70 to 2.8 microns 

if 

Halfwldth 

I 

ll nm I 

40 nm 
17 nm 
16 nm 


Center of Band <- 
380 nm 
500 nm 
875 nm 
940 nm 


TABLE III KADIOMETER CALIBRATION S: 
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TABLE V SOLAUTO PROGRAM (DDK) FEATURES 

USER 

DEFINED NAME AND FUNCTION OF KEY 
KEY # 

1 "Note to Tape" - Enter a user's note to be stored on tape. 

2 "Read ROMs" - Read all active ROMs, prompts for tape storage. 

3 "Photograph" - Activates camera, prompts for tape storage. 

4 "D to A" - Controls Digital to Analog feature of ROMs. 

5 "Intervals" - Sets photograph and ROM reading intervals. 

6 "Reset" - Resets all ROM and user defined functions. 

7 "ROMs On/Off" - Controls readings of ROMs two through five. 

8 "Time" - Displays time and prompts for reset. 

9 "Output Control" - Selects if records are to be written to tape and/or screen. 

10 "Photocount" - Lists current photo count and prompts for reset (36 maximum) . 

11 "EOT" - Locates end of data on tape and leaves tape positioned for writing. 

12 "A to D" - Controls Analog to Digital feature of ROMs. 

13 "Status" - Prints out list of current parameters in the form listed below: 
ROM^, i«2-5, 1 ~ enabled; 0 - disabled 

RDM DEL - Interrupt interval Cmlnutes) for reading ROMs 

PIC DEL - Interrupt Interval (minutes) for photographs 

ROM/TAPE - 1 - write records' to tape; 0 suppress 

ROM/DIS - 1 — write records to screen; 0 ~ not displayed 

PIC/TAPE - 1 — write records to tape; 0 — suppress 

Last ROM - Last ROM record 

Last PHOTO - Last photograph record 

PHOTO CNT - current photograph count 

TAPE REG - number of records written to tape 

Time - current time 

DELAY - current delay time to next Interrupt 

FLAG^ - values of flags: 1 - ROM record; 2 - Photo record; 3 - User inessage; 

4 - System -messages 

14 Pin 5 CONTROL - Output is high, low, low-^hlgh— low, high -low^hlgh, minimum 

pulse width Is 0.01 sec. 

15 "TAPE DUMP" - Reads tape by date and/or flag. 

16 "NEW TAPE" - Runs sequence for installing new data tape. 

17 "GRAPH" — Plots ROM data differentials vs. time by date. 

18 - Reserved for future use 

19 "Continue" - Resumes data acquisition after stopping program with UDK 20. 

20 "Stop" - Stops program execution until resumed by the use of UDK 19. 


TABLE VI DATA RECOVERY RECORD: HOURLY INTEGRATED GLOBAL SOLAR 
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Continuous global insolation measurements began at 1200 EST February 16, 1981. 
The strip chart record was too variable to allow recovery of these values. 
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* Continuous atmospheric emittance measurements began 00 EST February 24, 1981 
** Ice coated the silicon hemisphere o£ the sensor, 
i* The battery had to be changed. 


TABLE VIII DATA RECOVERY RECORD: HOURLY INTEGRATED DIFFUSE SOLAR 


m 

I 


I 

O 

m 


li 




o 

(O 


9^ 

W 


oR 

O w 
0) 

pii 0 ) 

u 

td (d 

a w 

(d 

P 


^'S 

S 8 

^ U 0 
4 ^ cn P 

0 (d 

H 

Sad 

1 


5 

d 




w 

43 




rs cr* 

) 1 

1 

1 

1 

1 

1 

1 


QO 

m 


rs o 

1 






O 

SO 

so 

•s 

so 









I I I • 0\ « 


I I I I I 


i I i t 


I I 


r -4 t> 
• • 
in m 


o 

o so 


IS 

CO 


§ ^ 


a 


0) 

sS 

M 



1 

ij 

0) 

*§ 

1 1 

p 

4J 

> 

a 

0) 

o 

o 

0) 

to 

o 

5Z5 

P 


m I I *4* <n 
o N 


I I I I I 


5 r 

<s 


I I I I I 


IS 

<1* 

« 

• 

OS 



■ 

OS 

• 

in 

• 

in 

Ol 

00 

in 

o 

o 

OS 

so 

so 

os 

Osi 

OS 

00 

o 

H 

o 

H 

OS 

OS 

so 


r> o 

tn 

M 

so 

M 

cr» 

cs 


so 

H 

OS 

H 


H 

rH 

CO 

CO 


•4 

H 

OI 

IN 

H 





so 

so 

r> 

f> 

IN 



& 

i 


cd 

B 

I 

Pm 


>d 

a 

I 


-& a 


0) 

I 


4 rH 

M 

i 00 

00 

) OS 

OS 

4 r-4 

iH 


8 

»a 


4J 

01 

tn 

1 

a 

3. 

4 

M 

^ j 

01 d 

■s ° 

5 1 
8 E 
^ S 

I 

u 

S 


I 

CO 

S 

nd 

M 

O 

3 

N 

4 J 

U 

'fJ 


o 

(0 

0) 

d 

o 


a* M 

•H O 

^ I 
0) 0) 


4 

+- 4 


31 




TABLE IX DATA BECOVERY RECORD: HOURLY mTEGRATED DIRECT SOLAR IRRADIANCE 
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TABLE X DATA BECOVERY RECORD EOR AI3TOMATED DATA ACQDISITON SYSTEM 

ONE AND TEN MINUTE INTEGRATED IRRADIANCES 
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TABLE XI AM, REGRESSION COEFFICIENTS FOR CLEAR SKT GLOBAL SOLAR IRRADIANCE 
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44. Normalized Cloudy Sky Global Solar Irradlance versus Cloud Cover Fraction 
for April 1982 

45. Normalized Cloudy Sky Global Solar Irradlance versus Cloud Cover Fraction 
for May 19.82 

46. Normalized Cloudy Sky Global Solar Irradlance versus Cloud Cover Fraction 
for June 1982 

47. Clear Sky Global Solar Irradlance versus Solar Zenith Angle for March 1981 
through February 1982 

48. Normalized Cloudy Sky Global Solar Irradlance versus Cloud Cover Fraction 
for March 1981 through February 1982 

49. Normalized Global Solar Irradlance versus Whole-Sky Fhotographlc Fraction 
for March 1981 through February 1982 
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50. Normalized Global Solar Irradlance vetaus Satelllte-PhotoprlnC Derived 
Fractions for Ylarch 1981 through February 1982 

51. Frequency of Occurrence of Specific Cloud Fractions obtained from the 
2926 LAFB Visual observations for March 1981 through Frebxuary 1982 

52. Frequency of Occurrence of Specific Cloud Fractions obtained from the 
95 Whole-Sky Photographic Values for March 1981 through February 1982 

53. Frequency of Occurrence of Specific Cloud fractions obtained from the 
53 Satellite Photoprint Values for March 1981 through February 1982 

54. Hampton Institute Study Region and Local Meteorological Data Sources 
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Figure 1. Automated Data Acquisition System Flow Chart 
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JULY 1981 



Figure 3 . Average Diurnal Variation of Global Solar Irradlance and 
Atmospheric (IR) Emlttance for July 1981. 
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Flgtire 16. Monthly Average Aerosol Extinction Coefficient, or Optical 

Depth, at 380 nanometers for the period Februai^ 1981 through 
June 1982. 
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Figure 17. Monthly Average Aerosol Extinction Coefficient, or Optical 

Depth, at 500 nanometers for the period February 1981 through 
June 1982. 




OPTICAL DEPTH 


♦ r 


gmOMAL^MK W 

SfwSquwjty 

Tt*'' ■ 


t 


57 



'' i ^ 

J i 


: mH OF YEAR 


Fl'gure»1'8 Extine tlen .Coefficient, or Optlcel 

Depth^^i «e.B75 nanoaeteri for the period fehruery 1981 through 
June 1982* 



ANGSTROM TURBIDITY COERFICIENT^ 30 


OMQINM. PAQE W 
OF POOR qUMilTY 


58 


iINGSTROH TURBIDITY COEFFICIENT, HONTHLY AVERAGE, FEB.8I-JUN.82 



NORTH OF YEAR 

Figure 19. Monthly Average Angstrom Turbidity Coefficient, B(9, for 
the measurement period February 1981 through June 1982. 
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Figure 54 Hampton Institute Study Region and Local Meteorological Pita Sources 
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DEVELOPMENT OF A SOLAR ENERGY MEASUREMENT LABORATORY FOR THE STUDY OF INSOLATION 
VARTATIONS AT HAMPTON, VIRGINIA. T. J. Griffin* . D. A. Whitney, and D. D. 
VentHbilc. Dept, of Physics and Engineering Studies, Hampton Institute^ Hampton, 
Virginia 23668. 

The purpoae of this three~year atudy is to investigate the cloud dependence 
of Incident solar rndiation at Hampton, Virginia, Solar irradiance at the 
Earth's surface la rolated to the extraterrestrial aolar irradiance, to radi- 
ation abaorbed cind emitted by the atmosphetc and clouds, and to radiation 
reflected by the Earth-atmosphere system. A ground-based measurement station 
h«e been establielied at Hampton Institute to monitor solar radiation, atmospher- 
ic emitted radiation, local cloud cover, and atsKispherlc turbidity. Continuous 
measurements of global, direct and diffuse solar radiation, and atmospheric 
Infrared radiation are made and stor«^d by computer, NOAA GOES-EAST satellite 
data are used to obtain albedo and .;ioud cover information. 

Interim analyses performed on the data include monthly averages of global 
insolation, infrared radiation, and atmospheric turbidity. Global insolation 
has been correlated with fractional cloud cover from March 1, 1981 through Feb- 
ruary 1, 1982 UBing the ARL empirical model. 

(Supported by NASA grant No. NAG- 1-87) 
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CUMPUTEH AUTOMATION OF A SOUAR RADIATION 
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A uu I Hr luidl ut l on maaBuramauti laboratory that inclMdaa laa 
l^raaiaiun apaotra pyranomatarii i praolalon Inirared radioMafcart 
a norimil Incidanaa pyrhtiiomttar and an all aky oamara haa been 
n\Hoiuai.«d no allow dlreon coinputer riedutiiitlon oC Inaolatlorv duna. 
Signula from the oplar Inetruinente are integrated and dliplayed on 
five dial n light -emit tlni‘*diode diapltiya, The integrated aignala are 
tranMferrod via a general porpoee interfaoe oard to u inierocomputerv? 
Wo have designed, contttrMCted and implemented hardware and Boftware 
ooiifigurationii that permit data aoduiaitioni Btoragei transfer i and 
display • System reiiabilltyteita have been performed and mean time 
between failure and ayatem down time have been character I sed, 
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